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Abstract. I describe recent work on inclusion of Pomeron loops in the high energy evolution. In 
particular I show that the complete eikonal high energy evolution kernel must be selfdual. 



Last year has seen renewed attempts to understand Pomeron loop contributions to 
the high energy evolution of hadronic cross sections in QCD. In recent years the study 
of the high energy scattering has centered around the so called JIMWLK evolution 
equation HI 0. Hfl . It describes the approach of the scattering amplitude to saturation 
due to multiple scattering corrections on dense hadronic targets, or in the diagrammatic 
language, the fan diagrams. The JIMWLK equation however only partially takes into 
account the processes whereby the gluons emitted in the projectile wave function at 
an early stage of the evolution, are "bleached" by subsequently emitted gluons, or the 
so-called Pomeron loops Jl, Hi, 0,01,01. 

Recently we have calculated corrections to the JIMWLK equation, which take into 
account some finite density effects in the projectile wave function (or equivalently, 
resum certain corrections away form the dense limit of the the target)[9]. We have also 
derived the evolution equation valid for dilute target, which is the opposite limit to that 
considered in JIMWLK||7j]. The most striking feature of the two results, is that they 
appear to be dual to each other. The improved JIMWLK equation is given by||9t] 
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where & denotes path ordering with respect to x~ and the field bf satisfies the "classical" 
equation of motion[9]. The low density limit evolution kernel (KLWMIJ) including the 
same type of corrections but in the target wave function derived in [7] is 1 : 

z klwmu += _ ^J{bf(zA[p]mzA[p])+b%zMPW(zMP]) 



1 Following the work 0], the same expression has also been obtained in 1 10] using the effective action 
techniques. 
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The two kernels are strikingly similar which suggests an intriguing duality between 
the high and the low density limits of the evolution kernel. 

In this contribution I follow [8] and show that indeed the full eikonal kernel for the 
high energy evolution must satisfy the property of self duality. The requirement that the 
evolution of the projectile and the target wave functions has the same functional form 
coupled with the requirement of Lorentz invariance of the scattering matrix, leads to the 
condition that the kernel of the evolution x[p, j^] must satisfy 
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where p is the charge density in the target wave function and a is defined byH! 

a a (x,x-)T a = ± 1 {x-y){u\y,x-)p a {y,x-)T a U{y,x-)), 



U(x,x~) = &>Qxp{i f dy~T a a a (x,y-)}. 
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I note that from the functional integral point of view this duality has been discussed 
earlier in lUlll . 

Consider the general expression for the S-matrix of a projectile with the wave function 
|P) scattering on a target with the wave function |T)||8|], where the total rapidity of the 
process is Y. The projectile is assumed to be moving to the left with rapidity Y — Yq (and 
thus has sizeable color charge density p~), while the target is moving to the right with 
rapidity Yq (and has large p + ). We assume that the projectile and the target contain only 
partons with large kr and k + momenta respectively: kr > A and k + > A + . The eikonal 
expression for the 5- matrix reads 

y Y = j Dp +a (x,x-)Wl[p + (x-,x)]Z p Y _ Yo [a], (5) 

where Z p is the S-matrix averaged over the projectile wave function 

Z p [a] = (P\ 6p e 'J dx ~ ! d2x P~ a ( x ) aa ( x y x ~) |p) . (6) 

where W T [a] is the weight function representing the target, which is related to the target 
wave function in the following way: for an arbitrary operator 0[p + ] 

(T| 0[p+(x)} \T) = J Dp +a W T [p + (x-,x)} 0[p + (x,x-)} . (7) 

The field a(x) is the A + component of the vector potential in the light cone gauge 
A = 0. This is the natural gauge from the point of view of partonic interpretation of the 
projectile wave function. In the formulae above we use hats to denote quantum operators. 



Note that the quantum operators p~ a (x) and p +a (x) do not depend on longitudinal 
coordinates, but only on transverse coordinates x. The "classical" variables a and p + 
on the other hand do depend on the longitudinal coordinate x~. This dependence, as 
discussed in detail in ED arises due to the need to take correctly into account the 
proper ordering of noncommuting quantum operators. Thus the ordering of the quantum 
operators /3 + in the expansion of O in the lhs of eq.0 translates into the same ordering 
with respect to the longitudinal coordinate x~ of p + (x~) in the expansion of 0[p + (x~)] 
in the rhs of eq.Q. 

As shown in [7] the functional W r [a] cannot in general be interpreted as probabil- 
ity density, as it contains a complex factor. This factor - the Wess-Zumino term, en- 
sures correct commutators between the quantum operators p a . In the present derivation 
we do not require an explicit form of this term, but the following property which is 
implicit in eq. © is crucial to our discussion. The "correlators" of the charge density 
(p ai (xi,x~[)...p an (x n ,x n )) do not depend on the values of the longitudinal coordinates 
xj , but only on their ordering 07J]. 

Note that one can define an analog of W T for the wave function of the projectile via 

(P\ 6\P~{x)\ \P) = J Dp a W p [p-} 0[p-(x,x-)] . (8) 

With this definition it is straightforward to see that E p and W p are related through a 
functional Fourier transform. To represent E as a functional integral with weight W p 
we have to order the factors of the charge density p in the expansion of eq.@, and 
then endow the charge density p~(x) with an additional coordinate t to turn it into a 
classical variable. This task is made easy by the fact that the ordering of p in eq.® 
follows automatically the ordering of the coordinate x~ in the path ordered exponential. 
Since the correlators of p (jc, ti) with the weight W p depend only on the ordering of the 
coordinates t\ and not their values, we can simply set t = x~. Once we have turned the 
quantum operators p into the classical variables p(x~), the path ordering plays no role 
anymore, and we thus have 

Z P (a) = J Dp a W P [p) gifdx- fd 2 xp a {x,x-)a"(x,x-)_ ^ 

We now turn to the discussion of the evolution. The evolution to higher energy can 
be achieved by boosting either the projectile or the target. The resulting S-matrix should 
be the same. This is required by the Lorentz invariance of the S-matrix. Consider first 
boosting the projectile by a small rapidity 8Y. This transformation leads to the change 
of the projectile S-matrix E of the form 

Substituting eq.dTOt into eq.© we have 




Where the second equality follows by integration by parts. We now impose the require- 
ment that the S-matrix does not depend on Yq 11211 . Since E in eq.© depends on the 
difference of rapidities, requiring that dy/dYo = we find that W should satisfy 

i!y wT = x[a '~L ]wT[p+] (i2) 

Thus we have determined the evolution of the target eq. (fT2l by boosting the projectile 
and requiring Lorentz invariance of the S-matrix. On the other hand the extra energy 
due to boost can be deposited in the target rather than in the projectile. How does W T 
change under boost of the target wave function? To answer this question we consider the 
relation between E and W together with the evolution of E. Referring to eqs.© and (flOb 
it is obvious that multiplication of if by a is equivalent to acting on W p by the operator 
—i8/8p, and acting on E p by 8/ 8 a is equivalent to multiplying W p by ip. Additionally, 
the action of ip and —i8/8p on W p must be in the reverse order to the action of 8 /8a 
and a on E p . This means that the evolution of the functional W p is given by 

^W p = x[-i^ip}W P [p}. (13) 

Although eq.([T3l refers to the weight functional representing the projectile wave func- 
tion, clearly the functional form of the evolution must be the same for W T . Comparing 
eq. (fT2l and eq.(fT3l we find that the high energy evolution kernel must, as advertised, 
satisfy the selfduality relation eq.©. This is the main result. 

The selfduality of the kernel is somewhat similar (although different in detail) to the 
duality symmetry of a harmonic oscillator Hamiltonian p — > x, x^ —p. One thus hopes 
that it may eventually be of help in solving the complete evolution equation, once it is 
derived. 
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